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ABSTRACT 



The UV and optical spectra from the HST FOS and from the FAST spectrograph 
at the Whipple Observatory reported by Kenyon et al (2001) at different phases are 
analysed, leading to new results about the configuration of AG Peg. The FAST spectra 
contain both Ha and H/3 lines, whose ratio changes with phase indicating that different 
nebulae contribute to each spectrum. In particular, the spectrum emitted from one of 
^vq \ the nebulae must be collision dominated, in order to justify the relatively high Ha/H/3. 

• Consistent modelling by the code SUMA, which accounts for the coupled effect of the 

photoionisation from the hot star and the shock, shows the important role of at least 
£NJ , three nebulae: 1) the nebula between the stars, downstream of the shock created by 

■ collision of the winds, which propagates in reverse towards the WD, 2) the nebula 

downstream of the shock expanding in the outskirts of the giant atmosphere, and 3) 
O ■ the shock propagating outwards the binary system, that is not reached by the hot 

source radiation. The relative contribution of the three nebulae to the final spectra 
at different phases are calculated for all the lines. It is found that the contribution to 
the different lines changes with the viewing angle of the system. The relatively low 
preshock density and magnetic field adopted in the modelling of the expanding shocks 
reveal that they are merging with the ISM, thus explaining the slow decline of AG 
Peg. 

' Key words: shock waves;stars: binaries: symbiotic-stars:individual: AG Peg 



1 INTRODUCTION 

AG Pegasi (AG Peg) is a symbiotic binary consisting of a 
white dwarf (WD) and a red giant (spectral type M3). A con- 
tinuing slow decline of the far-UV and near-UV continuum is 
noticed by Kenyon, Proga, & Keyes (2001). They claim that 
the eruption of the hot component in AG Peg may be the 
slowest classical nova outburst ever recorded. AG Peg evo- 
lution was followed by the observations (Kenyon et al 1993, 
Kenyon et al. 2001) and modeled by the interpretation of 
the line and continuum spectra (e.g. Contini 1997). 

The collision of the winds from the hot star and the red 
giant in symbiotic binaries creates a complex hydrodynamic 
structure (Nussbaumer 2000 and references therein). The 
emitting gas within the system is then ionised and heated 
both by the photoionising flux from the hot star and by 
shocks. In a previous paper (Contini 1997, hereafter Paper 
I) the observed continuum and line spectra from AG Peg 
at phases earlier than <j> = 7.12, have been explained by 
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composite models accounting both for photoionisation and 
shocks. Indeed, the complex structure of the system requires 
schematic models which represent roughly the real picture. 
The consistency of calculation results for the line spectra in 
the UV and in the optical ranges, as well as for the spectral 
energy distribution (SED) of the continuum leads to a better 
understanding of the physical characteristics of the system. 

Following the results of Paper I, two shock fronts are 
created at collision of the winds between the stars, one prop- 
agating in reverse toward the WD and the other propagating 
outwards the system. Moreover, recent observations of sym- 
biotic novae show that spectra are also emitted from down- 
stream of shocks propagating outside the system and even- 
tually merging with the ISM (Contini & Formiggini 2001). 
At phase 7.05 of AG Peg the expanding shock has already 
passed over the red giant and slightly accelerates propagat- 
ing in the atmosphere of the red giant, opposite to the WD, 
where the density gradient is decreasing. 

New observations are presented by Kenyon et al. (2001, 
hereafter KPK01). They provide optical and UV spectra at 
phases 7.34 - 9.44, acquired with Faint Object Spectrograph 



2 M. Contvni 



_ 1 1 1 1 1 1 1 1 1 1 1 


i i 1 i i i 


1 1 _ 




Ha 


\ - 




H/S 




Ha/H/3 




-/ - 




Hell 




\: 
V. 






T : 






\ - 


Hel 




\ - 




H7/H/3 










\ / 

1 1 1 1 1 1 1 1 1 1 1 


[FeVII] 

■ • ■ ■ 


i . 


8.4 8.6 8.8 


9 


9.2 



phase 



Figure 1. The evolution of some significant line intensities with 
phase. The lines are Hel 4471, Hell 4686, [FeVII] 6087, H/3, and 
Ha. 

(FOS) on board of the Hubble Space Telescope (HST) be- 
tween late 1993 and late 1996. Low resolution optical spectra 
with the FAST spectrograph mounted at the F.I. Whipple 
Observatory 1.5-m telescope are also reported. After an ac- 
curate analysis of the data, KPK01 conclude that "if the hot 
component emits enough high energy photons to account for 
the forbidden line fluxes, mechanical heating from the col- 
liding winds may explain the high electron temperature. A 
detailed photoionisation calculation which includes shock ex- 
citation would test this proposal." Responding to KPK01, I 
present further modelling of AG Peg on the basis of their 
observational data at the latest phases (<f> = 7.39 - 9.44), 
adopting for the calculations the code SUMA (Viegas & 
Contini 1994, Paper I). SUMA calculates the emitted spec- 
trum from the gas, accounting for the coupled effect of shock 
and photoionisation from an external source. The spectra in 
the different phases are considered in Sect. 2. The models 
and the modelling process are presented in Sects. 3 and 4, 
respectively. Results are discussed in Sect. 5 and concluding 
remarks appear Sect. 6. 



2 THE SPECTRA AT DIFFERENT PHASES 

Line fluxes observed by KPK01 were corrected for E(B- 
V)=0.1. There are several values for E(B-V) in the range 
0.08-0.15 obtained by different methods (Tomov, Tomova, 
& Raikova 1998, Penston & Allen 1985, etc.). Some signif- 
icant line fluxes and line ratios corrected for E(B-V) = 0.1 
are given in Fig. 1 as function of the phase. Fig. 1 shows that 
Hq/H/3 and H7 /H/3 at phases cf> < 8.8 are higher than case 
B values. A higher reddening correction of the spectra is not 
predicted because it would lead to very high values of the 
far UV line fluxes. Moreover, KPK01 notice that the depth 
of the 2200 A absorption feature does not vary with phase 



and is consistent with E(B-V) = 0.1. Finally, the best fit of 
the calculated continuum SED to the data of the IR dust 
bump at early phases (Paper I) corresponds to a relatively 
low gas-to-dust ratio, d/g = 2. 10~ 15 by number. 

We would like to investigate the origin of such high val- 
ues of Hq/H/3 and H7/H/?. We will refer to radiative trans- 
fer effects of HI, to collision excitation of H lines, and, more 
specifically, to the effect of the viewing angle of the emitting 
nebulae. 

Recall that densities in the emitting nebulae of symbi- 
otic systems are generally high (> 10 4 cm -3 ), sufficiently 
high that some self absorption occurs in the Balmer lines. 
This leads to a marked strengthening of the Ha line relative 
to the rest of the Balmer series. In case B, as the optical 
thickness of Lya increases, the first effect in the decrement 
occurs when H/3 converts into Ha (Cox & Mathews 1969). 
Particularly, for large optical thickness of the Lya line, the 
increase of the Ha optical thickness leads to an increase of 
both Ha/H/3 and H7/H/3 (Cox & Mathews 1969, Fig. 1). 
In fact, Ha is mostly scattered because any L/3 photons it 
forms are quickly absorbed and converted back to Ha, while 
H/3 is absorbed and converts to Ha+ Pa. 

Moreover, Drake & Ulrich (1980) claim that Ha/H/3 ra- 
tios greater than 10 can be obtained with optical depths of 
the order of 100 or more by Balmer self-absorption alone. 
Particularly, conditions which tend to produce very steep 
Balmer decrements are found to be large optical depths in 
Ha, low photoionisation rates, and values of the electron 
density between 10 s cm~ 3 and 10 12 cm -3 (see also Krolik 
& McKee 1978). 

In the case of symbiotic systems, collisional processes 
are very effective, because shock fronts are created by col- 
lision of the winds between the stars with velocities ^ 40 
kms -1 . Such velocities correspond to temperatures > 20,000 
K in the postshock region. Higher velocities are actually ob- 
served leading to higher temperatures. Collisional excitation 
increases with temperature and prevails for T ^ 20,000 K. 
Collision strengths are oc ETj 2 exp(- Ey/kT) (Contini & Al- 
drovandi 1983) where E^ is the energy of the level. Radia- 
tion from the hot star, on the other hand, cannot heat the 
gas to ^ 20,000 K. So, high Ha/H/3 line ratios are predicted 
in the spectra emitted from shock dominated nebulae. 

It was shown in Paper I that the spectra of AG Peg 
at different epochs are emitted mainly by two nebulae : one 
in the downstream region of the shock front propagating in 
reverse towards the hot source and the other downstream of 
the shock front expanding outwards in the atmosphere of the 
red giant opposite to the hot source. In particular, the gas 
is collisionally heated in the outer region of the expanding 
nebula, while is heated by radiation in the internal region 
of the same nebula facing the hot source. So, collision dom- 
inated spectra will prevail at photometric minimum. Fig. 1 
shows that the Ha/H/3 line ratio has a rough maximum at 
phase ~ 8.5, which could be due to collisional effects. We 
will come back to this point after modelling the spectra. 

When the system evolves through the photometric min- 
imum and the hot source is occulted, the Hell lines that are 
emitted by a photoionised gas, should strongly decline. On 
the contrary, Fig. 1 shows that the Hell 4686 flux does not 
change with phase more than about 22% throughout the 
period between phases 8.33 and 9.44. This behavior reveals 
that nebulae heated and ionised by radiation contribute to 
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the spectra near maximum as well as near photometric min- 
imum. 



3 THE MODELS 

SUMA calculates the emitted spectrum from a gas in a plane 
parallel geometry. The nebula between the stars is ionised 
and heated by both black body radiation from the WD and 
the (reverse) shock. Therefore, radiation reaches the very 
front edge of the shocked nebula. On the other hand, radi- 
ation and shock act on opposite edges of the nebula down- 
stream of the expanding shock which propagates in the red 
giant atmosphere. 

The input parameters are : the shock velocity, V s , the 
preshock density, no, the preshock magnetic field, Bo, the 
colour temperature of the hot star, T* , the ionisation param- 
eter, U, the geometrical thickness of the emitting gaseous 
nebula, D, and the abundance of He, C, N, O, No, Mg, Si, 
S, A, and Fe relative to H. For all models T, = 10 5 K. 

We calculated a large grid of models with the following 
criteria: to fit the most significant line ratios and to be con- 
sistent with the modelling reported in Paper I which was for 
earlier photometric phases. 

Three models were selected from the grid, ml, m2, and 
m3. Model ml represents the nebula in the downstream re- 
gion of the shock propagating in reverse toward the WD. 
The downstream nebula, facing the red giant, is however 
close to the WD and the photoinising flux from the hot 
source reaches the shock front edge of the nebula. Model 
m2 represents the nebula downstream of the shock prop- 
agating outwards the system in the red giant atmosphere 
opposite to the WD. The downstream nebula toward the 
WD is ionised and heated by both the shock and radia- 
tion from the hot source. Model m3 shows the same shock 
characteristics of model m2, but the emitting nebula down- 
stream is not reached by the WD radiation flux. Model m2 
is radiation dominated (RD), while m3 is shock dominated 
(SD). Model m3 being SD, the emitted spectrum is mainly 
excited by collision. Recall that at later phases the shock 
front, merging with the ISM, is disrupted by R-T and K- 
H instabilities (see Contini & Formiggini 2001). Therefore, 
the ionising flux from the hot source will be eventually pre- 
vented from reaching the nebula in the downstream region 
of the shock by intervening matter. 

So, in the following, the "reverse shock" will refer to 
model ml, the "RD expanding shock" to model m2, and 
the "SD expanding shock" to model m3. Models ml and m2 
are matter-bound, while m3 is radiation bound (see Sect. 
4.2). 



4 MODELLING THE SPECTRA 

The spectra contain lines in a relatively large range of ion- 
isation levels and in the different wavelength domains. Re- 
call that the distribution of the fractional abundances of the 
ions downstream (Fig. 2) depends on the input parameters 
and that the line intensities increase with the volume of gas 
emitting the corresponding ion and with the square of the 
density. 

The choice of the models depends on line ratios from 



different ionisation levels and is constrained by the general 
model of the AG Peg system. 

Between cf> = 7.05 and <f> = 8.07 the conditions of the 
reverse shock (between the two stars) could not change sen- 
sibly because the reverse shock is most likely a standing 
shock, so, we adopt input parameters at <j> = 8.07, V s =60 
kms _1 and no = 10 6 cm~ 3 similar to those found at = 7.05, 
V s =40 kms^and n = 3 10 6 cm" 3 ( Paper I, Table 4, col- 
umn 2). In particular, Bo = 10~ 3 gauss is adopted for ml, 
which is suitable for isolated giants (Bohigas et al 1989). 

The preshock density of the shock propagating through- 
out the atmosphere of the red giant in the present mod- 
elling (10 3 cm~ 3 ) is lower than that calculated at phases 
7.05 - 7.12 (Paper I, Table 4) by a large factor (~ 3000). 
The shock velocity is slightly lower (100 kms _1 at 0=8.07 
and 150 kms _1 at = 7.05). A lower n indicates that the 
shock has proceeded towards the outskirts of the giant at- 
mosphere. In about 820 days (between =7.05 and =8.07) 
the shock expanding with a velocity of ~ 100 kms _1 could 
reach a distance of ~ 7. 10 14 cm and is propagating outwards 
the system, merging with the ISM. Therefore, the preshock 
density of model m2 is low and the magnetic field is also 
reduced to 10 -4 gauss. Both no and Bo are, nevertheless, 
still higher than found in the ISM (1-10 cm~ 3 , and 10~ 5 - 
10~ 6 gauss, respectively). 

The relative abundances adopted at ^ 8.07 are C/H 
= 3.3 10~ 4 , N/H = 5.3 10~ 4 , O/H = 9.6 10~ 4 , Ne/H = 8.3 
10" 5 , Mg/H = 2. 10~ 5 , Si/H = 3.3 10~ 5 , S/H = 1.6 10~ 5 , 
A/H = 6.3 10~ 6 , and Fe/H = 4. 10~ 5 , as were found at 
0=7.05-7.12 (Paper I, Table 4). Mg/H is, however, similar 
to that at = 2.74 (Paper I, Table 1) because at = 7.05 
Mg lines were not observed and the Mg/H abundance ratio 
could not be determined by modelling the corresponding 
lines. The adopted relative abundances show that N/H and 
O/H ratios are higher than solar. Relatively high values of 
N/H and O/H are characteristic of symbiotic systems. 

The modelling process of the spectra is roughly de- 
scribed in the following. Recall that model ml represents 
the reverse shock, m2 the RD expanding shock , and m3 
the SD expanding shock. The results corresponding to sin- 
gle models are summed up adopting relative weights at each 
phase and the final results are compared with the observa- 
tions. The relative weights which lead to the best fit of the 
observed spectra in the different epochs are related to the 
configuration of the system and are discussed in the next 
section. 

High ionisation-level lines (e.g. NV, [MgV], [FeVII], etc) 
increase with the ionisation parameter U and with V s . In 
fact, the temperature downstream behind the shock front 
is oc V s . NV/CIV ratios calculated by model ml with rel- 
atively low velocities (V s < 80 kms -1 ) are lower than ob- 
served, therefore, they must be counter-balanced by line 
ratios higher than observed, calculated by the m2 model. 
The NV/CIV ratio calculated by SD models (m3) with V s ~ 
100 kms _1 is very low. So, the only way to obtain a high 
NV/CIV is to adopt a high U in model m2. This leads to 
high Hell/CIV line ratios. 

The NV/CIV ratios calculated by a RD model with 
ionisation and shock acting on opposite edges, are not high 
enough, if the density is high (10 6 cm -3 ). Also, at high 
densities, the ratio of the intermediate ionisation-level lines 
(NIII/CIV, Silll/CIV, and CIII]/CIV) would be higher than 
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Tabic 1 

Line ratios to CIV 1550 from selected models 



line 
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m2 
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0.15 
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0.1 


1.0 
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HI 4340 


0.04 
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observed by a factor > 10. High densities are thus rejected 
for model m2. This is consistent with the previous discussion 
(Sect. 3), namely, that the RD expanding shock propagates 
in a low density medium. 

Low ionisation-level line ratios (e.g. CII] /CIV and 
Mgll/CIV) are low for models ml and m2 (Table 1). Low 
ionisation-level lines are high in SD models, therefore, model 
m3 must be accounted for. The relatively high CII] /CIV and 
Mgll/CIV ratios calculated by m3 counter-balance in the 
averaged spectra the low values found by models ml and 
m2. 

A large grid of models is calculated in order to achieve 
the best fit of all the line ratios. 

4.1 Calculated line ratios 

The line ratios calculated by models ml, m2, and m3 are 
presented in Table 1. The input parameters are given in 
the bottom of the table. The parameter s, which appears 
in the last row of Table 1 distinguishes between the two 
cases: radiation from the star reaches the shocked edge of 
the emitting nebula (s=0) and radiation and shock act on 
opposite edges (s=l). The calculated absolute H/3 flux and 
the CIV/H/3 line ratios are also given in order to calculate 



the flux of the lines. Recall, however, that the spectra are 
calculated at the nebula but are observed at Earth. The 
adjusting factor depends on the distance to Earth of AG 
Peg, as well as on the distances of the nebulae to the WD. 

4.2 Distribution of the different ions throughout 
the nebulae 

In order to better understand the line ratios calculated by 
the different models, the distribution of the fractional abun- 
dance of some interesting ions as well as of the electron den- 
sity and temperature downstreams is shown in Fig. 2 for 
models ml (top diagram), m2 (middle diagram), and m3 
(bottom diagram). The shock front is on the left; black body 
radiation reaches the left edge in the top diagram (ml) and 
the right edge in the middle diagram (m2). The middle di- 
agram is divided in two equal symmetric parts in order to 
give equal importance to the shock dominated region (left) 
and to the radiation dominated region (right). Model m3 
(bottom diagram) is shock dominated (U=0.). In all the di- 
agrams the ions are represented by the following symbols. 
Solid lines : H +1 ; dotted lines : He +2 ; short-dashed lines: 
N +5 ; dash-dotted lines : C +4 ; long-dashed lines : C +3 ; long- 
dash-dotted lines : C +2 . 

Fig. 2 shows that the temperature throughout the neb- 
ulae never decreases below 10 4 K downstream of the reverse 
shock (top diagram). This is due to the coupled effect of a 
relatively high U and a relatively small geometrical thick- 
ness of the nebula, too small to host a large recombination 
region. 

The profile of the temperature is different in the RD ex- 
panding nebula (middle diagram). The gas behind the shock 
front is heated to T > 1.5 10 5 K by the shock correspond- 
ing to V s = 100 kms -1 . The high ionisation level lines (e.g. 
NV) are strong due to the coupled effects of the relatively 
high temperature and the strong radiation flux from the hot 
source. The cooling rate (oc n 2 ) of the gas is sped up by free- 
free, free-bound, and line emission, therefore a large cool re- 
gion appears in the left side of the diagram. Notice that H 
remains fully ionised and the N +5 /N fractional abundance 
is high in the shocked side of the nebula, even if the tem- 
perature drops, because they are maintained by the strong 
radiation flux from the hot source. C +4 prevails throughout 
the nebula corresponding to the ml and m2 models owing 
to the high U. 

The shock dominated model (m3) shows different condi- 
tions, characteristic of SD models. A relative large region is 
dominated by the C +2 ion. The absolute flux of the H/3 line 
is low because the emitting region throughout the nebula is 
relatively small. 

Interestingly, Ha/H/3 is higher in the postshock region 
(left side of the middle diagram) of model m2 than in the 
region reached by the radiation from the hot source (right 
side) and is definitively high throughout all the nebula rep- 
resented by model m3. 



5 RESULTS 

The best modelling of AG Peg should, indeed, imply a huge 
number of models in order to represent the complex nature 
of symbiotic systems (see Nussbaumer 2000). The models, 
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Figure 3. Top panel : A schematic representation of AG Peg 
system. Bottom panel : The orientation of the binary system axis 
at different phases. 



however, are schematized to three main types representing 
the three emitting nebulae described in the previous sec- 
tions; therefore, the results show only a consistent approach. 



5.1 Comparison of calculated with observed line 
ratios 

The calculated spectra at each phase are compared with 
observations in Table 2. Each spectrum results from the 
weighted sum of the three models. The relative weights of 
the models are shown in the last rows of Table 2. The weights 
are very different for different models because they depend 
on the emitting area of the corresponding nebula. 
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Table 2 

The comparison of calculated and observed line ratios to CIV 1550 
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01 


0.05 


0. 


.21 


0.11 


[NcV] 3426 


0. 


.022 


0.13 





.073 


0.35 





.023 


0.30 





.15 


0.74 


HI 4340 





.056 


0.04 





.07 


0.06 





048 


0.058 





.113 


0.110 


Hel 4471 





.004 


0.0001 





007 


0.0001 





.005 


0.0001 





.010 


0.0002 


Hell 4686 





.06 


0.08 





.159 


0.14 





.131 


0.129 





.273 


0.250 


[FeVII] 6087 1 






0.06 





.019 


0.15 





012 


0.128 








w(ml) 






1.7(-3) 






1.3(-4) 






2.(-4) 






0.0 


w(m2) 






20. 






20. 






20. 






20. 


w(m3) 






6.67 






3.0 






3.0 






2. 



1 [FeVII] flux from the FAST spectra at 0=8.36 and 0=9.27 



A very schematic representation of the system and the 
angles of view at the different phases are given in Fig. 3 (top 
panel and bottom panel, respectively). 

At — 8.07 we see almost all of the shock fronts, inter- 
nal and external to the system (Fig. 3) . The distance of the 
reverse shock from the white dwarf, ri , is of the order of the 
binary separation (in Paper I it is found that the front is 
very close to the red giant surface), so, we adopt ri =5 10 13 
cm ; the radius of the RD expanding front is 12 (0=8.07)= 
r2(0=7.O5)+A r, where A r is the distance covered by the 
shock front between phases 8.07 and 7.05. Then, X2 = 3.2 
10 15 cm (Paper I) + 7 10 14 cm (Sect. 4) is ~ 4 10 15 cm. 

Actually, (r 2 /ri) 2 ~ 6400 at = 8.07 is lower than 
w(m2)/w(ml) = 1.2 10 4 by a factor < 2. The ratios of the 
weights which appear in Table 2 are determined by the best 
fit of the calculated to the observed data. The discrepancy 
by a factor of < 2 indicates that the radiation emitted by the 
reverse shock is slightly absorbed by the outer atmosphere 
of the red giant. At later phases, = 8.33 and = 9.27, 
the radiation intersects the red giant at angles that depend 
on the phase angle, and the path length throughout the red 
giant, 1 = 2 trq cos a (trg is the radius of the red giant 
and a is the angle between the binary system axis and the 
viewing direction from Earth), is higher at 0=8.33 than at 
0=9.27. So, at 0=9.27 the radiation from the reverse shock 
is less absorbed than at 0=8.33, in agreement with the trend 
of the w(m2)/w(ml) ratio that is 10 5 at 0=9.27 and is 1.5 



10 5 at 0=8.33 . The trend is clearly depending on the phase 
angle, however, the trend does not depend only on the angle 
because the radius of the expanding shock is not symmetric 
around the system (see Nussbaumer 2000, Fig. 1). 

Regarding the SD shock, the ratio of the relative weights 
w(m3)/w(m2) does not follow any specific trend, because it 
depends on fragmentation of matter and on the filling factor. 
Also those ratios were chosen by the best fit to the data. 

Considering that the AG Peg system is very complex, 
we have focused on the most significant lines, e.g. NV 1240, 
CIV 1550, Hell 1640, CIII] 1908, CII] 2325, etc. in the UV, 
and Hell 4686 in the optical range. NIV] 1486 is overesti- 
mated by the calculations. OI 1305, Hel 4471 and Mgll are 
underestimated by a large factor. The contribution from a 
model corresponding to a large low temperature region could 
improve this result. [NeV] 3426/ CIV line ratio is overpre- 
dicted. Notice, however, that this ratio refers to lines emit- 
ted from one or two different nebulae, indicating that more 
model components should be considered. The same is valid 
for the [FeVII] 6087 /CIV line ratio. In this case, however, 
some iron could be locked into dust grains. 

A slightly lower N/H could improve the fit of the NV 
1240 and NIV] 1486 line ratios to CIV, particularly at phase 
9.44. N/H approaching to the ISM value, and the relatively 
low no and Bo adopted for models m2 and m3, confirm that 
the merging process is taking place. 

In the averaged spectra we found Hq/H/3< 3, lower than 
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the observed line ratios, because the weights of ml and m2 
models are relatively high. A possible explaination is given 
by Shore (1992) who claims that in a moving medium with 
an internal velocity gradient, the optical depth of the line At 
~ Av/(dv/dl). Av=Avd, where Avd is the Doppler width 
of a line with rest frequency vo- The reduction of the veloc- 
ity gradient causes a rapid increase in the line opacity , even 
if the density has dropped enough as to make the equiva- 
lent static line profile completely optically thin. Indeed, a 
velocity gradient is predicted in the outer atmosphere of 
the red giant. Adopting V s = 150 kms _1 at <f> — 7.12 (Con- 
tini 1997), V s =100 kms _1 at the latest epochs, and dl=10 14 
cm, At should increase by a factor > 10 7 , leading to higher 
Ha/H/3 (see Sect. 2). 

We can now compare the luminosities of the lines cal- 
culated by the models with those observed by KPK01. Con- 
sidering, for instance [NeV] 3426 at <f> = 8.07, the calcu- 
lated luminosity is ~ 3 10 30 erg s _1 while the observed one 
(KPK01, Table 4) is 4 10 31 erg s _1 adopting a distance to 
Earth, d=800 pc. This line is mainly emitted from the neb- 
ula corresponding to model ml. A better fit to the observed 
luminosity could be obtained adopting no=3 10 6 cm~ 3 in- 
stead of 10 6 cm -3 for ml. The preshock density used at <f> 
= 7.05 (Paper I, Table 4) was 3 10 6 cm" 3 . 

The distance of the SD expanding nebula from the cen- 
ter of the binary system, r3, at <^>=9.44 is obtained compar- 
ing the luminosities of calculated and observed CII] lines. 
From CII] o6s x d 2 = CII] m3 x r|, r 3 results ~ 10 16 cm. The 
CII] line is considered because it is emitted mainly from the 
outer SD nebula at photometric minimum. The distance is 
larger by a factor > 2 than the distance of the RD expand- 
ing nebula corresponding to model m2. Even if the RD and 
SD expanding shocks are at different distances from the sys- 
tem center, they may belong to the same wave. In fact, the 
"fingers" created by R-T and K-H instabilities at the shock 
front can be elongated up to a large fraction of the shock 
front radius (see Contini & Formiggini 2001). 

In summary, the interpretation of the data by a rel- 
atively small number of models is sound, even if approxi- 
mated. Anyway, some interesting results follow. 

5.2 Relative contribution of the different nebulae 
to the spectra 

The relative contribution of the emitting nebulae to each 
line is shown in Table 3. The weight of models m2 and m3, 
representing the nebulae propagating outwards in the atmo- 
sphere of the red giant, relative to the weight of the neb- 
ula corresponding to the reverse shock between the stars 
increases with the angle determined by the phase fraction 
(.07, .27, .33, and .44), indicating that the external nebula 
is not present around the whole system, but only in the re- 
gion opposite to the WD. 

The nebula downstream of the reverse shock dominates 
at </> = 8.07 except for a few lines (Silll], CIII], CII], and 
Mgll). A large contribution to these lines comes from the 
SD expanding shock. The RD expanding shock contributes 
partly to NIV, [NeV], and to [FeVII]. 

The relative contribution of the three nebulae to the 
lines are similar at <j) = 8.33 and at <j> — 9.27, the reverse 
shock being less dominant and indicating that the relative 
contributions of the nebulae depend on the viewing angle of 



the system rather than on the phase. At both phases 8.33 
and 9.27 the RD expanding shock is dominating for most 
of the lines (NV, Hell 1640, NIV, [MgV], [NeV], Hell 3203, 
[NeV] 3426, HI 4340, Hell 4686, [FeVII], and Ha). The SD 
expanding shock dominates the low ionisation level lines. 

At <f> = 9.44 the contribution from the reverse shock 
becomes negligible because the nebula between the stars is 
occulted by the red giant. The high ionisation level lines 
are produced mainly by the RD expanding shock and low 
ionisation level lines come from the SD nebula. 

The relative contribution of the nebula between the 
stars to the Hell 1640 line decreases monotonically with 
phase fractions .07, .27, .33, and .44 On the other hand, the 
relative contributions of m2 to NV 1240 are 14.8%, 61.4%, 
68.75%, and 99.0%, at the same angles, respectively. 

Concluding, the evolution of the spectra is more likely 
related to the viewing angle of the system, not particularly 
to the phase between <f> = 8.07 and (j> = 9.44. We do not 
have data relative to viewing angles higher than 0.44, which 
could hopefully strengthen the results. 

5.3 Evolution of the line fluxes with phase 

It is now possible to understand the trend of the lines shown 
in Fig. 1. Table 3 shows that the contribution from the SD 
shocked nebula (m3) to Ha is only a few % , therefore, the 
calculated Ha/H/3 line ratio is likely to be ~ 3 throughout 
the whole period between phase 8.07 and 9.44. On the other 
hand, the observed value is > 5 towards photometric mini- 
mum. This high value may be explained considering that the 
RD expanding shock is represented by a composite model 
which accounts for both the photoionisation by the WD and 
the shock. The side of the emitting nebula facing the ob- 
server at photometric minimum is shock dominated (Fig. 2, 
middle diagram). So, in this phase we see mainly Ha and 
H/3 fluxes from this side of the nebula, which are produced 
by collisional ionisation. Moreover, the optical depth in the 
radiation dominated region is high (see Sect. 4.1) leading to 
high Ha/H/3. 

Hell 4686 and [FeVII] do not show a periodic change 
with phase because a large fraction of these lines is emit- 
ted by the nebula downstream of the RD expanding shock. 
So, the decrease of the reverse shock contribution at phases 
approaching the photometric minimum is compensated by 
the large contribution from the RD expanding shock. These 
lines show, therefore, less impressive variability up to <j> ^ 
9.2. At 4> = 9.27 the calculated F A 4686 is 0.11 ergcm~ 2 s _1 , 
while at <j> — 9.44 Fa4686 = 0.06 erg cm -2 s _1 . The decrease 
of the calculated flux between these two phases is in agree- 
ment with the observed decrease shown in Fig. 1 after <j> = 
9.27. 

Hel 4471 is emitted mainly by the SD expanding shock 
(m3) and the flux is very low, Fa4471= 2 1 ~ 5 ergcm~ 2 s _1 . 
Therefore, this line flux shows a maximum in Fig. 1 when 
the system is near maximum. The same is valid for Ha and 
H/3. Notice, however that Hel 4471/CIV line ratio is under- 
predicted by the models. 

All the lines show an asymmetric behaviour before and 
after maximum. This is explained by the fact that the emit- 
ting nebulae are not symmetric relative to the axis joining 
the center of the WD with the center of the red giant. In the 
hydrodynamical picture of symbiotic systems presented by 
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Table 3 

Relative contributions (in %) to the line intensities 







0=8.07 






0=8.33 






0=9.27 








0=9.44 




line 


ml 


m2 


m3 


ml 


m2 


m3 


ml 


m2 


m3 


ml 


m2 


m3 


NV 1240 


84.7 


14.78 


0.54 


30.11 


68.75 


1.13 


37.60 


61.38 


1.02 








99.0 


0.99 


SilV] 1394 


78.90 


3.68 


17.40 


34.4 


20.98 


44.61 


44.66 


17.7 


37.64 








43.69 


56.31 


SilV.OIV] 1403 


91.35 


5.39 


3.25 


50.45 


38.96 


10.59 


61.04 


30.64 


8.32 








85.86 


14.14 


NIV] 1486 


92.9 


5.25 


1.81 


53.94 


39.88 


6.19 


64.3 


30.9 


4.79 








91.41 


8.59 


CIV 1550 


92.27 


3.57 


4.16 


56.47 


28.56 


14.97 


66.62 


21.9 


11.48 








75.89 


24.11 


Hell 1640 


86.94 


13.05 


0.014 


33.74 


66.23 


0.032 


43.92 


56.05 


0.026 





.0 


99.97 


0.029 


OIII] 1664 


67.78 


1.35 


30.87 


25.39 


6.60 


68.01 


34.36 


5.81 


59.83 


u. 


n 

u 


13.81 


86.19 


NIV 1719 


79.63 


20.36 


0.003 


23.02 


76.98 


0.0051 


31.51 


68.49 


0.0045 








100. 


0.0 


NIII] 1750 


64.84 


1.35 


33.8 


23.04 


6.29 


70.66 


31.54 


5.60 


62.86 


0. 


.0 


12.81 


87.19 


Silll] 1892 


30.21 


0.64 


69.14 


6.78 


1.89 


91.32 


10.0 


1.82 


88.11 








3.30 


96.70 


CIII] 1908 


47.39 


0.81 


51.74 


13.07 


2.93 


84.0 


18.78 


2.74 


78.48 








5.44 


94.55 


CII] 2325 


0.27 


0.002 


99.73 


0.05 


0.004 


99.95 


0.071 


0.004 


99.93 





.0 


0.0074 


99.99 


[MgV] 2783 


81.05 


18.94 


0.0035 


24.65 


75.34 


0.006 


33.48 


66.50 


0.0056 








99.99 


0.0052 


Mgll 2793 


4.09 


0.033 


94.88 


0.90 


0.077 


99.02 


1.38 


0.077 


98.54 








0.128 


99.87 


[NeV] 2976 


85.60 


14.37 


0.032 


31.28 


68.65 


0.069 


41.19 


58.75 


0.058 








99.94 


0.06 


Hell 3203 


87.0 


12.98 


0.013 


33.88 


66.09 


0.031 


44.08 


55.89 


0.026 








99.97 


0.028 


[NeV] 3426 


73.3 


26.66 


0.0043 


17.38 


82.61 


0.006 


24.46 


75.54 


0.0056 





.0 


100. 


0.0 


HI 4340 


86.7 


12.93 


0.37 


33.61 


65.56 


0.836 


43.78 


55.5 


0.71 








99.23 


0.77 


Hel 4471 


69.80 


2.99 


27.2 


25.95 


14.56 


59.49 


35.03 


12.77 


52.2 


0. 





28.76 


71.24 


Hell 4686 


85.94 


14.09 


0.009 


31.78 


68.20 


0.012 


41.75 


58.23 


0.017 





.0 


99.98 


0.0177 


H/3 4861 


100. 


0.000 


0.0 


100.0 


0.0 


0.0 


100.0 


0.0 


0.0 








99.99 


0.007 


[FeVII] 6087 


72.84 


27.16 


0.0015 


17.02 


82.98 


0.002 


23.98 


76.0 


0.002 





.0 


100. 


0.0 


Ha 6563 


85.77 


12.79 


1.43 


32.8 


63.98 


3.22 


42.89 


54.38 


2.74 


0. 





97.04 


2.96 



Nussbaumor (2000, Fig. 1 left) the wake is, in fact, strongly 
asymmetric, so that the emitting nebulae could partly disap- 
pear in some location at angles corresponding to high phase 
fractions. 



6 CONCLUDING REMARKS 

The analysis of the spectra from the HST and from the 
FAST spectrograph at the Whipple Observatory reported 
by Kenyon et al (2001) at different phases leads to new re- 
sults about the configuration of AG Peg. The FAST spectra 
contain both Ha and H/3 lines, providing the tool for mod- 
elling the system at different phases. Interestingly, the trend 
of the periodical variability of Ha/H/3, as well as the trend 
of the most significant lines indicate that different nebu- 
lae contribute to the final spectra. In particular, the high 
Ha/H/3 line ratios (~ 5) suggest that shock dominated 
models must be accounted for. Moreover, the trend of the 
Hell 4686 line variability with phase indicates that radiation 
dominated nebulae are strongly contributing to the spectra 
also near the photometric minimum. 

Consistent modelling by the code SUMA, which ac- 
counts for the coupled effect of photoionisation from the hot 
star and shocks, shows that at least three nebulae provide 
the spectral components : 1) the nebula between the stars, 
downstream of the shock created by head-on collision of the 
winds, which propagates in reverse towards the WD, 2) the 
nebula downstream of the head-on-back shock expanding in 
the outskirts of the red giant atmosphere and reached by the 
photoionising flux from the WD, and 3) the nebula down- 



stream of the same shock front which is shock dominated, 
i.e. it is not reached by the photoionising radiation flux. In 
fact, the matter in this region is very clumpy because the 
shock front is disrupted by R-T instability, and radiation 
from the WD is prevented from reaching the region by the 
intervening clumps in some locations. 

The fit of calculated to observed line ratios is acceptable 
enough to justify the choice of the models, so that the rela- 
tive contribution of the different nebulae to the final spectra 
at different phases can be calculated for all the lines. It is 
found that the contribution of the nebula between the stars 
prevails for most of the lines near maximum (0 = 8.07), 
declining towards higher phases and even disappearing at 
phase 9.44. Low ionisation level lines (e.g. CII] and Mgll) 
are mainly emitted from the SD expanding nebula. The evo- 
lution of the spectra between phases 7.34 and 9.44 is more 
likely related with the viewing angle of the system, and not 
to the temporal evolution of the system. 

Models m2 and m3 which represent the expanding shock 
are characterized by a relatively low preshock density and 
magnetic field, much lower than found by modelling AG Peg 
at ^ 7.12 (Paper I). This suggests that the shock, expand- 
ing in the red giant atmosphere opposite to the WD, has 
already reached the outskirts of the system and is merging 
with the ISM. At <f> = 9.44 the expanding shock has swept up 
a large portion of matter with relatively low density. The line 
and continuum fluxes emitted from downstream the shock 
are oc n 2 and, therefore, decrease as more IS matter is swept 
up, leading to the slow decline of AG Peg continuum and 
line fluxes observed by KPK01. 
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